Preprint typeset in JHEP style - PAPER VERSION 



Precise evaluation of the Cut Constructible part of 
one loop amplitudes within D-dimensional unitarity. 



Achilleas Lazopoulos 

Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland 
E-mail: lazopoulosOitp . phys . ethz . ch 



Abstract: A method to separate pentagon contributions from the evaluation of the cut 
constructible part of primitive amplitudes within the framework of D-dimensional unitarity 
is proposed. The cut constructible part is thus reconstructed with significantly higher 
accuracy as demonstrated by the numerical analysis presented here. 

Keywords: NLO computations, QCD. 



1. Introduction 



Precise predictions for signals and backgrounds at hadron colliders are the necessary in- 
gredient for the unveiling of new physics. Such predictions can only be achieved with 
calculations beyond the tree-level leading order in the perturbation expansion over the 
strong coupling constant. 

Next to leading order QCD calculations consist of a virtual part and a real emission 
part. The real emission is nowadays automated either in the Catani-Seymour dipole ap- 
proach [1], implemented in [2], [3], [4, 5], [6], [7], or the Frixione-Kunszt-Signer approach[8, 9] 
recently implemented in [10]. The virtual part that involves one loop integrals has been 
tackled in many different ways, most of which involve some kind of reduction to master 
integrals. Each one loop integral is expressed as a linear combination of a small number of 
simpler integrals which are known in a closed analytic form. Traditionally this reduction 
is performed analytically for each one-loop Feynman diagram appearing in the calcula- 
tion and the coefficients are known as functions of the kinematical variables of the process. 
There are many generic frameworks implemented that have been used the last years in very 
demanding calculations (see for example [11, 12, 13] used in [14, 15, 16] or the publically 
available [17] used in [18, 19]). 

Alternative approaches based on unitarity have been employed extensively in the past 
to recover one loop amplitudes (see e.g. [20, 21]). Following a generalization of the unitarity 
idea to quadruple cuts [22] and a novel approach in reduction that allows one to partially 
fraction numerically one loop integrals at the integrand level on a per point base [23, 24], the 
numerical evaluation of one loop amplitudes using tree-level quantities as building blocks 
appeared to be possible and efficient enough to tackle processes with many final state 
partons. Since then there have been mainly three approaches developed along the lines 
of numerical reduction with tree-level objects as building blocks: the Black Hat approach 
[25, 26](based also on the reduction method of [27], implemeneted in [26, 28]), the D- 
dimensional unitarity approach [29, 30, 31, 32] (implemented in various independent codes 
[33, 34, 35, 36, 37] ) and the helac-lloop approach[24, 38, 39, 40, 41, 42, 43] (implemented 
in [44] ) . All three approaches have already delivered differential cross sections for processes 
of prime phenomenological importance [45, 46, 47, 48, 49, 50, 51, 35, 52, 53, 54]. 

Due to their infrared and ultraviolet singular nature, loop amplitudes need to be evalu- 
ated in D = 4 — 2e dimensions, and also the degrees of freedom of unobserved particles must 
be set according to some regularization scheme. Within the framework of D-dimensional 
unitarity, the virtual amplitude is considered as a function of the dimensionality of the 
virtual particles circulating in the loop. Keeping the dimensionality integer, one can eval- 
uate both the coefficients of the cut constructible (CC) part and the full rational part of 
the one loop amplitude employing D-dimensional cuts that decompose the coefficients into 
tree-level amplitudes in D-dimensions. The price paid for the simplicity of the algorithm 
is the necessary inclusion of pentuple-cuts (absent in other approaches). 

Pentuple-cut residues involve the loop momentum in a high power and can pose a 
threat to the numerical accuracy of the result when evaluating the reduction coefficients for 
phase space points of degenerate kinematics (when any Gram determinant of the external 
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momenta approaches zero). The percentage of such points varies with the process and 
the kinematical cuts on the external partons employed to avoid soft and collinear singular 
regions. As a rule of thumb, their number increases in processes with many gluons and few 
fermion pairs. 

In particular, pentuple-cut residues when large, are seen to reduce the numerical ac- 
curacy of the rational as well as the cut-constructible part of the amplitude. On the other 
hand, it is known [29] that the CC part can be fully reconstructed by restricting the loop 
momentum (and therefore the whole unitarity method) in four dimensions. It is the aim 
of this work to show how the pentuple-cut residues can be decoupled from the evaluation 
of the CC part which is then evaluated with much higher precision. 

2. D-dimensional unitarity 

We will follow the version of D-dimensional unitarity presented in [31] and the notation of 
[37]. We employ color decomposition to write the virtual amplitude as a linear combination 
of gauge-invariant primitive amplitudes [21]. 

Each such primitve amplitude is a sum of one-loop integrals that need to be dimen- 
sionally regularized. Apart from the dimensionality of the integration, that is continued 
from L> = 4toD = 4 — 2e there is also the dimensionality, D s of the internal, unobservable, 
particles, which depends on the particular regularization scheme: it is set to 4 in the Four 
Dimensional Helicity scheme (FDH) [55] or to 4 — 2e in the 't Hooft - Veltman scheme 
(HV) [56]. The dependence of the amplitude on D s is linear (with the exception of fermion 
loops). Seen as a function of D s , 



and one can evaluate numerically A a for two values of D s = Di,2, determine the D s inde- 
pendent Aq, A\ from A Dl and A D2 , and thus get the full A Ds for arbitrary D s . Thereafter, 
setting D s to 4 recovers the FDH scheme and to 4 — 2e the 't Hooft - Veltman scheme. 
During all this the dimensionality of the loop momentum is kept to arbitrary D with the 
constraint D < D s . Only after the reduction is performed one can set D — ► 4 — 2e and 
evaluate the master integrals as a series in e. The problem of numerically evaluating full 
one loop amplitudes (including the rational part) reduces then to the problem of evaluating 
A Dl ' 2 . In order to accommodate fermions, one has to choose D\ 2 = 6, 8. 

Each such primitive amplitude A Dl ' 2 is partially fractioned in terms of a complete set 
of master integrals. 



where = {i\,i2, ■ ■ ■ ,ik} a se t of k indices corresponding to k of the propagators that 
appear in the primitive amplitude (</& runs over the entire set of possible subsets of k 



Aq + A x ■ D s 



(2.1) 




(2.2) 
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propagators among those that appear in the primitive) . We parametrize the unintegrated 
coefficients e(Z) , d(l) , c(l) ,b(l) as follows 

ej 5 (i) = ej 5i0 

dj 4 (l) = dj 4fi + dj^iai + dj 4 , 2 (4 + ^4,301^5 + dj 4 ^a\ 

cj 3 (0 = cj 3 ,o + cj 3j iai + cj 3j2 a 2 + cj 3i3 (ai - a\) + aia 2 (cj 3i4 + cj 3 , 5 ai + cj 3j6 a 2 ) + 

2 1 2 , 4 

cj 3 ,7aia 5 + cj 3j8 a2a 5 + cj 3i9 a 5 

&J 2 (0 = & J 2 ,0 + &J 2 ,l a l + &J 2 ,202 + 6/2,303 + 0j 2 ,4(a? - al) + &/ 2 ,5(a 2 ~ a 3) + h J2fi a l a 2 + 

bj 2 jaia^ + bj 2:8 a 2 a 3 + bj^gOg 

(2.3) 

with a,i = l-rii and n« are the unit vectors that span the subspace that is orthogonal to the 
physical subspace spanned by the momenta involved in the propagators 

Once all coefficients are recovered the primitive amplitude can be written as 

A = A cc + A R (2.4) 
A cc = dj Afi Ij A + Y cj 3 ,o!j 3 + E b J2,o J J2 ( 2 - 5 ) 

Ja J3 J2 

where dj 4: o includes dj 4 $ as well as contributions from pentagon scalar integrals that are 
reduced to boxes. 

J& J3 J2 

3. System of equations and Pentagon contamination 

In order to evaluate the unknown coefficients one has to solve an OPP-like system: begining 
from the pentagon, for every a set of specific loop momenta lj k]S for which Dk(lj k - S ) = 
for all k e Jk is determined. For those loop momenta, the corresponding propagators arc 
then on-shell, and the pentuple-cut coefficient e(l s ) can be determined as the residue of the 
loop integral, which is seen to be equal to a product of tree-level color-ordered amplitudes. 
Given the pentagon coefficient, one can proceed to the boxes. The box coefficients d(l) is 
equal to the residue minus a counter term from pentagons. Then one proceeds to the triple 
cuts and so on and so forth. This is what is meant by the statement that the system of 
OPP equations is triangular (see [23, 24]). 

The system of equations we have to solve is, therefore, 



E C h,r9r(ls) 
r=0...9 



e h,s = = e h,0 

M U D ' h/h{ls) -0...4 

c h,s - m h - 2^ n . . : run ~ 2^ 
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dj 4 (U 




cj 3 (U 



bi 2 ,rh r {l s ) 



UkeJi/h D k(h) 



Dj 3 /i 2 (ls) 



r=0...9 



(3-1) 



where is the value of the residue of the loop amplitude with respect to the propagators 
indicated by Ik and the choice of loop momentum l s . Moreover, 



The right hand side of the above equations is just a shorthand for the parametrization of 
eq.2.3, i.e. the functions f r (l),9r(l), h r (l) are just the quantities multiplying the coefficients 
di 4 ,r,ci 3 , r ,bi 2 ,r in eq.2.3. 

Among the unknown coefficients {eo, d r , c r , b r } the subset {eZo,i, Co.. .6, &o...s} can be 
evaluated separately, by choosing the loop momentum l s in four dimensions. In that case 
{^3,4,5^7^,9,69} drop out. The subset {do ; ±, Co.. .6, &o...s} reconstruct the four-dimensional 
part of the amplitude and <io,co,6o are the coefficients of the master integrals in the cut 
constructible part of the amplitude. Therefore their value should be independent of the 
value of pentuple cuts. This however is not explicit in the above system, where pentuple 
coefficients are involved as counterterms in the equations for quadruple-, triple- and double- 
cut coefficients. 

That this is an undesirable situation is made clear from studies of numerically unstable 
points. Pentuple-cut coefficients are carrying the effect of possible small Gram determinants 
to the fifth power. In the framework of the Vermaseren- van Neerven base [29], [57], Gram 
determinants appear as denominators in the evaluation of the loop momentum part that is 
fixed to solve the unitarity constraints of the cut. Phase space points for which the Gram 
determinants are approaching zero lead to numerically large components of the pentuple- 
cut loop momentum and therefore to large pentagon coefficients. Their effect must cancel, 
however, when solving the four-dimensional subsystem of eq. 3.1 since the cut-constructible 
master integral coefficients can be derived by the purely four-dimensional of [29] where 
there are no pentagons. 

The cancellation of the pentagon contribution within the cut-constructible sub-system 
is the main source of numerical instabilities in the calculation of the CC part of the ampli- 
tude. One can see the correlation of numerically unstable phase space points (as detected 
by the inaccurate reconstruction of their pole coefficients in the e = 2 — D/2 expansion 
of the amplitude) with the appearance of small Gram determinants in the pentuple cut 
or equivalently with the appearance of large pentagon coefficients. This has already been 
noted in [34]. 

It is also shown in the plots of fig.l for the case of six gluons. Very similar correlation 
patterns are found for more external particles and for amplitudes with external fermions. 

The effect of pentagon contamination in the calculation of the CC part of the amplitude 
is exhibited in fig. 2-3 for the six gluon primitive and the qq — ► 4g leading order primitive. 
The percentage of points that can be classified as numerically unstable based on the pole 
accuracy criterion varies from 5% downwards, assuming reasonable final state cuts (e.g. a 



J a \i b = Ja ■ h — h 
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Figure 1: Correlation of the maximum pentagon coefficient with the double and single pole accu- 



racies defined in the text e 2 ,i = log 10 
the phase space pointset. 
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See section 5 for details about the generation of 
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Figure 2: Double and Single pole accu- 
racy in ./V = 6 gluons with the canonical 
approach. See section 5 for details about the 
generation of the phase space pointset. 



Figure 3: Double and Single pole accuracy 
in qq — > 4g with the canonical approach. See 
section 5 for details about the generation of 
the phase space pointset. 



cone separation angle of R = 0.7 and a px cut of 20GeV). If one would like to explore the 
collinear or soft limits, of course, the number of unstable points would increase. 

4. Pentagon isolation 

One can adopt the attitude of completely discarding the problematic points (assuming 
that they don't populate particularly contingent regions of the phase space) or resorting 
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Figure 4: Double and Single pole in N = Figure 5: Double and Single pole in N = 
7 gluons with the canonical approach. See 8 gluons with the canonical approach. See 
section 5 for details about the generation of section 5 for details about the generation of 



the phase space pointset. 



the phase space pointset. 



to quadruple precision (with 32 or 64 digits in the floating point arithmetics) to evaluate 
them accurately. We would like to suggest in this paper a method to reduce the number 
of unstable points drastically. 

Given that the system 3.1 is linear, we split it in two systems of equations, isolating 
the counterterms that include pentagon coefficients. 
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(4.3) 

The b-subsystem contains exclusively contributions coming from the pentagon either di- 
rectly or through box and triangle counter terms, while the a-subsystem has contributions 
from quadruple-, triple- and double-cuts as well as counter terms coming from them. 

When the loop momenta are chosen in 4-d the b-subsystem coefficients of the boxes 
Dq 1 are non-zero. In that case there is no freedom in the choice of the box loop momenta. 

There are two solutions of the box unitarity constraints, lj 4 . ± = Vi A ± yJ—V^n^. The 

coefficient d b l4 is then 

4,=i<4<w + 4«-))=-i£«* Got + jot) (44) 

However all pentagon master integrals are reduced to boxes, so the final coefficient of each 
box master integral including contributions from the reduction of the pentagons is 

d h , 0J uii = d a hfl + d b hfi + d h:0 (4.5) 

and it is easy to see that 

4,„ - -<, - \ £ + (4.6, 

J 5\I 4 

This means that the box coefficient <ir 4) o,/«zz can be evaluated without the contribution of 
the pentagons when the loop momenta are in 4-d. Furthermore, we have 

i>j B // 4 (i-); 

(4.7) 

from which we get 

^('^--E^U^^ + wnJ (4 ' 8) 



with 

*j 8 (0 = -^= (4-9) 
The b-subsystem for the triple coefficients is then 



-b ej 5 ^Ja/iS 1 ^ 



^ 3 n fc6 j 6/ / 3 ^(^) ^ 3 D kl (i s ) 

J 5|/ 3 J 4|J 3 J 5|J 4 



(4.10) 



-7- 



The summation symbols can be interchanged since S,/^ = Ej 4 | Z Sjr| 7 > so we 

have 



4, 
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ej 5 
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n 



fceJ 5 // 3 
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J 4|/ 3 



1 f l + Zj 4 (Z s ) 1 



D JJL/Ia (l a )2 \D J5/J4 (l + ) D J5/JA {1_) 



(4.11) 



The quantity inside the brackets vanishes identically when l s is in 4D and it solves the 
unitarity constraints of the 1% triple cut. Note that this cancellation is a kinematical 
identity and is therefore independent of whether the external or internal particles are 
fermions or gluons. The cancellation has been numerically checked to hold up to 12 digits 
in the double precision version of the code for processes with up to 8 external particles of 
all flavors. 

This means that all c^ 3 s vanish identically if l s is in 4D, and so do the triangle coeffi- 
cients c b l3 6 for all 13 . Similarly, the bubble coefficients b b l2 8 vanish identically. 

Therefore the whole b-subsystem doesn't contribute to the CC part of the amplitude. 
It is, however, necessary for the evaluation of the rational part. 
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Figure 6: Double and Single pole in N = 6 
gluons with the new approach. See section 5 
for details about the generation of the phase 
space pointset. 



Figure 7: Double and Single pole in qq — > 
Ag with the new approach. See section 5 for 
details about the generation of the phase 
space pointset. 



In this setup, one can thus disentangle completely the pentagon coefficients from the cut 
constructible part of the primitive amplitude. The difference this makes in the numerical 
stability of the cut constructible part is rather spectacular as can be seen in figs. 6, 7, 8, 9, 
where the double and single pole accuracy plots of fig.2,3,4,5 are revisited. The CC part 
is now evaluated by solving the two subsystems separately. In figs 14-15 we show the 
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Figure 8: Double and Single pole in N = 7 
gluons with the new approach. See section 5 
for details about the generation of the phase 
space pointset. 



8 gluons ++--++-- 




Double and Single pole accuracy 



Figure 9: Double and Single pole in N = 8 
gluons with the new approach. See section 5 
for details about the generation of the phase 
space pointset. 



accuracy of the CC part of the finite result for a primitive with gluons and the leading 
color primitive with one fermion pair, as compared with the quadruple precision result. 

The improvement is significant and the percentage of points when one needs to go to 
higher precision in order to reconstruct the CC part are in the worst of cases at the per 
mille level. 

The importance of reconstructing the CC part of the amplitude with a numerically 
accurate method is amplified by the fact that, in the majority of phase space points, the 
CC part gives the larger contribution to the finite part of the amplitude. In fig. 10-11 we 
show scatterplots of the total finite result versus the ratio of CC/ finite for the gg — ► 4g 
primitive for two different helicity configurations, and in fig. 12-13 for the qq — * Ag leading 
color primitive. The percentage of points where the ratio deviates from 1 by, say, 10% 
depends on the particular helicity configuration, but it is clear that for the majority of 
points (from 75% for pure gluonic primitives to 95% for primitives that include fermions) 
the CC part is the larger contribution to the finite result for the primitive amplitude. 



5. The rational part 

When the loop momentum l s is five-dimensional, the coefficients d\ s ,c\ s ,b\ s are non-zero, 
and moreover depend on the numerically dangerous pentagon coefficients linearly. But 
knowledge of their explicit functional form allows us to reconstruct them in a way that 
partially avoids numerical cancelations. 
To this effect one notes that 
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Figure 10: Scatterplot of the absolute value 
of the total finite result for the primitive am- 
plitude, \FIN\ versus the ratio \CC\/\FIN\ 

for the + H h+ helicity configuration of 

the 6 gluon gluonic primitive. See section 5 
for details about the generation of the phase 
space pointset. 



Figure 11: Scatterplot of the absolute value 
of the total finite result for the primitive am- 
plitude, \FIN\ versus the ratio \CC\/\FIN\ 

for the + + H helicity configuration of 

the 6 gluon gluonic primitive. See section 5 
for details about the generation of the phase 
space pointset. 



and, since for the box we have Is = + a\n^ (for which I 2 = 0), 

(h + Qk) 2 = Q 2 + 2Z S ■ q k = q 2 + 2V ■ q k + a x n ■ q k = x + a x y (5.2) 



So 



from which 



$ =- y — ek n (5.3) 



keJ 5 /i 4 



4 4 ,0= E g * 2 2 h 2 ( 5 - 4 ) 

d h,i= J2 *fc 2 2 U , (5.5) 
keJ 5 /h /4 ^ 4 /4 

with etj 4 = J-VjI- Similar expressions can be found for ^2,3,4. The box counter term 
functions d b j 4 (l) (needed for the triple-cuts) can then be reconstructed for any I, with the 
pentagon coefficient pulled out as a common factor. This prevents cancelations of large 
terms while solving the system. 

However, in the majority of the problematic points, the rational part is inaccurately 
reconstructed due to the inaccurate evaluation of the the pentuple-cut residue. The above 
treatment therefore, only partially improves the numerical accuracy of the rational part, 
as shown in figs. 18-19 where the accuracy of the rational part relative to the rational 
part calculated with quadruple precision is shown with and without the pentagon isolation 
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Figure 12: Scatterplot of the absolute value 
of the total finite result for the primitive am- 
plitude, \FIN\ versus the ratio \CC\/\FIN\ 

for the + H 1 helicity configuration of 

the qq — > Ag leading color primitive. See sec- 
tion 5 for details about the generation of the 
phase space pointset. 



6 gluons ++--++ 




Figure 13: Scatterplot of the absolute value 
of the total finite result for the primitive am- 
plitude, \FIN\ versus the ratio \CC\/\FIN\ 

for the + H 1 h helicity configuration of 

the qq — > Ag leading color primitive. See sec- 
tion 5 for details about the generation of the 
phase space pointset. 



qqgggg ++--++ 




Figure 14: CC in the canonical and the new 
approach pole in N = 6 gluons. See section 5 
for details about the generation of the phase 
space pointset. 



Figure 15: CC in the canonical and the 
new approach pole in qq — > Ag.See section 5 
for details about the generation of the phase 
space pointset. 



approach described above. We expect that evaluating the pentuple cuts only, in higher 
precision, would remedy most of the problematic points left. 

Phase space points are selected uniformly by an implementation of the RAMBO[58] 
algorithm, with cuts E T > 10~ 2 ^/s,\n\ < 4,Ai? = s/An 2 + A<p 2 > 0.7. The definition of 
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8 gluons ++--++-- 





Figure 16: CC in the canonical and the new 
approach pole in N = 7 gluons. See section 5 
for details about the generation of the phase 
space pointset. 



Figure 17: CC in the canonical and the new 
approach pole in N — 8 gluons. See section 5 
for details about the generation of the phase 
space pointset. 



6 gluons ++--++ 



qqgggg ++--++ 





Figure 18: The total finite part with the 
canonical and the new approach in N = 6 
gluons. See section 5 for details about the 
generation of the phase space pointset. 



Figure 19: The total finite part with the 
canonical and the new approach in qq — > Ag 
gluons. See section 5 for details about the 
generation of the phase space pointset. 



accuracy that we adopt is the standard one: 

\x N -x P \ 

e = lo gio i i (5-6) 

\ X P\ 

where x can be the coefficient of the double or single pole in the e-expansion of the primitive, 
or the CC part or the total finite part, xjy is the numerical value of x and xp is the predicted 
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value, either from the analytic expression for the poles [59, 60, 61] or from numerical 
evaluation in quadruple (64 digits) precision, for the finite part. 

6. Computational efficiency improvement 

Thanks to the improved evaluation of the cut constructible part of the amplitude, the 
number of points for which one needs to resort to quadruple precision arithmetics is reduced 
drastically. In practice one can avoid quadruple precision in 20 — 35% of the previously 
problematic points. Since evoking quadruple precision typically involves a cpu penalty 
factor of the order of 100, the improvement can be significant in terms of total cpu time, 
particularly in the case of pure gluonic amplitudes or amplitudes with one fermion pair. 
For example, in the case of the qq — > gggg used as an example above, with the tight cuts 
we have adopted in this work, with the canonical approach one has 3.6% of the phase space 
points with relative accuracy of less than 10 -4 . The new approach brings the percentage 
down to 2.3%. This translates into a 30% reduction of the necessary cpu time needed to 
evaluate all points with relative accuracy better than 10~ 4 . 

7. Conclusions 

We have suggested a new method to disentangle the contribution of pentuple-cuts from the 
evaluation of the Cut Constructible part of one loop amplitudes within the framework of D- 
dimensional unitarity without increasing the computational complexity of the algorithm. 
We have showed that in all cases studied the Cut Constructible part is now evaluated 
with high precision. Knowledge about the explicit functional form of the contributions 
to the rational part that originate from pentuple-cuts allows for a small improvement in 
its evaluation. As a result, the percentage of phase-space points for which one needs 
to resort to higher precision arithmetics (which is computationaly expensive) in order to 
reconstruct the virtual amplitude correctly is drastically reduced. This leads to a significant 
improvement in the efficiency of the D-dimensional unitarity approach without the need 
for restricting the phase space available for final state jets or vetoing of problematic phase 
space points. 
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